Metastable gravity on classical defects 
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We discuss the realization of metastable gravity on classical defects in infinite-volume extra di- 
mensions. In dilatonic Einstein gravity, it is found that the existence of metastable gravity on the 
defect core requires violation of the Dominant Energy Condition for codimension N c — 2 defects. 
This is illustrated with a detailed analysis of a six-dimensional hyperstring minimally coupled to 
dilaton gravity. We present the general conditions under which a codimension N c > 2 defect admits 
metastable modes, and find that they differ from lower codimensional models in that, under certain 
conditions, they do not require violation of energy conditions to support quasi-localized gravity. 

PACS numbers: 04.50.+h, ll.10.Kk, 98.80.Cq 



I. INTRODUCTION 



In recent years, there has been great interest in 
models of infinite-volume extra dimensions. The 
Dvali-Gabadadze-Porrati (DGP) model realizes four- 
dimensional (AD) gravity on a membrane in a flat 
five-dimensional (5D) bulk by adding a large, induced 
Einstein-Hilbert term to the worldvolume of a three- 
brane The Einstein-Hilbert term on the brane is 
generically induced by quantum loops of matter fields on 
the brane, and realizations of this mechanism are known 
in more fundamental theories like string theory 0, 0, 13 • 
One of the main phenomenological attractions of DGP 
is its celebrated "self-accelerated" cosmological solution, 
giving small late-time acceleration without the need to 
invoke dark energy [{| . Higher codimensional realizations 
of DGP gravity are known to be non-trivial, because a 
straightforward generalization of the modelgives singular 



propagation and other difficulties 



ic model gives sing 



On the other hand, codimension N c > 1 theories 
have very attractive properties, foremost because of their 
unique way of addressing the cosmological constant prob- 
lem [H, Q • Being non-local in the infrared, these models 
evade Weinberg's no-go theorem, unlike compact extra 
dimensional theories. The idea put forward in for ex- 
ample is that, in infinite- volume extra dimensional theo- 
ries, one can have AD nearly flat solutions even though 
the brane tension is large ( "the AD cosmological constant 
curves the bulk, not the brane"). In more detail, using 
the induced gravity setting, it was argued that there are 
solutions of supercritical branes in infinite bulk with ten- 
sion T > Mf that support a regular geometry if one 
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allows for an inflating worldvolume with Hubble rate: 

, 4N l/(N c -2) 



H ~ AL 



Ml 
T 



(1) 



M* is the fundamental higher dimensional Planck mass 
(constrained to be quite low ~ 10~ 3 eV) and N c the num- 
ber of codimensions. In predicting an inflation rate in- 
versely proportional to the vacuum energy density these 
models provide an interesting way of solving the cosmo- 
logical constant problem. 

In Refs. 0, 0, 0, El different ways of realizing 
DGP-like gravity in higher dimensional context were dis- 
cussed. A fundamental difference with the codimension 
one model is that AD gravity is mediated there by massive 
resonant graviton states, with the width of the resonance 
much smaller than its mass ^1 (DGP allows AD gravity 
by a broad resonance peaked around zero mass). How- 
ever, the general task remains to find a physical frame- 
work where gravitational dynamics can be approximated 
by a DGP-like action 



S = Mi 



f2+7V c 



g\Rd N ' +4 X. (2) 



Here g and R are the determinant and scalar curvature 
of the induced metric on the brane while g and R are the 
same quantities for the bulk metric. 

In analogy with the action J2J), in Ref. [Tlj . 4Z?-like 
gravity is realized by a "regularized" action: 



S = M. 



2+7V c 



g\RF(X)d N ° +i X 



(3) 



where J 7 is a sharply peaked profile around the origin 
of the transverse space, which can be interpreted as a 
varying Planck mass in the transverse dimensions. This 
action has the same features as Eq. in that gravity 
in the brane core is weakly coupled while the bulk is 
strongly coupled. It was shown that by appropriately 
choosing the profile T , one propagates AD gravity by 



2 



resonances in the brane core 1 . Indeed, the Newtonian 
potential between localized sources on the brane, at a 
dimensionless distance x from eachother, can be obtained 
by considering the exchange of Kaluza-Klein modes 0, 

MM 



roc -Mi 

V s {x)ot / p(M) dM. 

Jo x 



(4) 



Here p(M) is the spectral density on the brane, which is 
modified from its empty space value in the presence of 
quasi-bound states. For instance, in flat 51? spacetime, 
one recovers a static potential varying as l/x 2 from a 
constant spectral density. One can consider the case in 
which the spectral density is, in addition to its standard 
behavior, strongly peaked around a particular mass M c 

p{M) = 1 + 5(M-M C ). (5) 

From Eq. > this would lead to a gravitational potential 



V g (x) oc — | c 



-M c x 



(6) 



and one recovers 4D-like gravity on the scales 1 -C x <C 
1/M C , which is physically attractive for M c <C 1. Inter- 
estingly, such systems predict generically an infrared as 
well as ultraviolet modification of gravity. In the regular- 
ization of |Tlj this behaviour was indeed recovered, using 
appropriate regulating profiles for the action Eq. JSJ. 

From Eq. lf2"|. it is clear that the previous mechanism 
could also be realized by letting y/\g\ vary transverse to 
the brane instead of the Planck mass. This has been 
studied in Refs. 0, 0] in the context of warped extra 
dimensions, with or without asymptotic flatness. 

Following the previous discussing, we want to focus 
in this paper on the following question [T^. Ecj : is it 
possible to realize AD gravity by (meta)stable states in 
the context of classical field theories? Topological de- 
fects serve as natural branes in non-perturbative field 
theory [2ll I22I I23I Wa\ EsT ] and they have been inten- 
sively studied in the context of warped extra dimen- 
sions [1^,|23,|2^,|2^|, while an obvious setting to gene rate 
a varying Planck mass is dilatonic gravity |30ll3lll3a . l33| . 
The formulation of these theories using an underlying 
non-linear sigma model realization is important to ad- 
dress their (supersymmetric) phenomenology 0] and in- 
ternal consistenc y ( for debate on the DGP model see 
e.g. Refs. 

Indeed, the action © suggests that in an appropri- 
ately formulated theory of dilatonic gravity, one might 
achieve metastable states in the core of a topological de- 
fect if the dilaton condenses strongly around the core and 



1 The varying gravitational coupling constant denned by T is anal- 
ogous to a varying dielectric constant in electrostatics. The ap- 
pearance of resonant states is supported by this analogy. 



falls off sharply outside the defect. Clearly, such a model 
would involve both variations in the metric through the 
defect energy-momentum tensor, and variations of the 
"bulk Planck mass" by means of the dilaton. It is not 
clear a "priori how these effects are balanced, and one 
has to construct explicit solutions to make quantitative 
statements about their properties. 

This paper is organized as follows. We first construct 
an illustrative example of a local hyperstring in six di- 
mensional dilatonic gravity. Considering a massive dila- 
ton, we find that the dilaton naturally condenses around 
the hyperstring core, in a region set by its Compton wave- 
length, and that the geometry of the system is asymp- 
totically flat. The full non-linear solution is computed 
numerically. Next we study the propagation of AD ten- 
sor perturbations on this background and find that there 
are no metastable states appearing in this setting. 

We then ask under which general conditions (on the 
defect forming matter and gravity content) tensor (with 
respect to the AD brane) resonant states exist in the de- 
fect core in infinite-volume extra dimensions. It will ap- 
pear that a violation of the Dominant Energy Condition 
(DEC) is a necessary condition for this to happen in codi- 
mension N c = 2. 

In the last section we restate the same question for the 
case of codimension N c > 2 defects. It is shown that, at 
least under reasonable conditions, metastable gravity is 
allowed without violation of positive energy conditions. 
We leave however explicit constructions for future work. 

Finally we conclude, and stress that our results 
only concern classical field theoretical realizations of 
metastable gravity. As mentioned before, there exist a 
number of different ways towards realizing metastable 
gravity on branes in infinite- volume extra dimensions, 
using the DGP loop induced effects or explicit construc- 
tions in string theory. Our work evidently does not ad- 
dress those models. 



II. 4D HYPERSTRING COUPLED TO 
DILATONIC GRAVITY 

In the following we consider the action for a four- 
dimensional hyperstring coupled to scalar-tensor gravity 
in a six-dimensional spacetime [30L l3lT ] : 



S= f ^[R-g^d^d^-Um^^x 



(7) 



Md b x, 



where g AB is the 6D metric with signature (— , +,+,...), 
R its Ricci scalar, <fi the dilaton field with a potential 
U (</)), = 327r 2 G 6 /3, G 6 being the 6-dimensional grav- 
ity constant. In order to allow for topological vortex 
configurations, we include a complex scalar field $: 

Anat = -\g AB CP A $) f T> B H> ~ - jH AB H AB , (8) 
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in which capital Latin indexes A,B... run from to 5 
and H AB is the electromagnetic-like tensor defined by 



Hab — d A C B d B C A , 



(9) 



where C B is the 1-form connection. The U(l) covariant 
derivative T> A is defined by 



V A = d A - iqC A 



(10) 



where q is the charge. The potential of the scalar field 
$ is chosen to break the underlying U(l) symmetry and 
thereby allow for hyperstring configurations, 



A 



V(<S>)='i{m 2 - V 2 f 



(11) 



A. Background Solution 

In terms of the dimcnsionlcss coordinate 

p = m^r = VXr/r, (17) 

mh being the mass of the Higgs boson, the Einstein tensor 
reads 



2/3.. 3. 2 il) 3.w 
= g^m h [-a + -a +- + -a- 



Go 



s / 3 . n „ . UJ 

ml ( -a 2 + 2a- 



u) 2 m\ [ 2ij + ^a' 



(18) 



where A is a coupling constant and r\ = (\&\) is the magni- 
tude of the scalar field vacuum expectation values (VEV) . 

The variations of the action (7J with respect to the 
metric and the dilaton field lead to the equations of mo- 
tion in the scalar-tensor gravity sector, in the Jordan 
frame one gets 



G AB = e-K 2 T AB + 28 A 



t;9ab [3<9 a 



V x d x c 



U(cj>)]+V A d B cj>-g AB V x d x 



(12) 
(13) 



where G AB is the TV-dimensional Einstein tensor and T A 
the stress tensor of the hyperstring: 



5g* 



9ab£-u 



(14) 



The metric is chosen so as to respect the cylindri- 
cal static symmetry in the two extra dimensions and 
Poincare invariance along the remaining four spacetime 
coordinates; 



where the dot stands for differentiation with respect to 
p. From Eqs. J3J) and (|14fl . the stress tensor reads 



/ 1 , 9 1 ip 2 Q 2 , , 1 m 2 Q 2 
T, v = -g, v ( ^mfo 2 + + V(<p) + 



T, 



™2 -2 



2 uj 2 

1 V 2 Q 2 , lm 2 Q 2 

2—- V[ip) 



Tee = u 



2 • 2 



VQ 2 



2 q 2 uj 2 ' 

ImlQ 2 



2 uj 



V{<p) 



2 q 2 uj 2 



In terms of dimensionless quantities, the Einstein and 
dilaton equations of motion 112H and l|13|) can be recast 
into 



—a H — a 
2 2 



w 3 . w 

1 — a — = ac 

ZD 2 ZD 



-r- 



2 PQ 2 



ZD~ 



ezu z 



3 



3 w 
-a H 

2 m 



1 U 
2ml 1 



(20) 



d.s 2 = g AB dx A dx B = e^rj^dx^dx" + dr 2 + uj 2 (r)d9 2 , 

(15) 

where r/^ is the four dimensional Minkowski metric of 
signature (—,+,+,+), and (r, 6) the polar coordinates 
in the extra dimensions. Greek indexes /i, v . . . run from 
to 3 and describe the brane worldvolume. 

We chose the Nielsen- Olesen ansatz for the matter field 
leading to a n hy perstring configuration with unit winding 
number [29ll38|. 



$ = ^(r)e* 9 , C e = -[1-Q(r)} 



<1 



(16) 



The only non-vanishing component of the electromag- 
netic tensor is He r = Q 1 1 q- The dilaton is also assumed 
to respect the cylindrical symmetry and depends only on 
the extra radial coordinates </>(r). 



w 





p f 2 Q 2 



















^(/ 2 -l) 2 + 

1 u 
"2< 



2<T + -a 2 = ac 



3 



f Q i 



1 u 



2 mi 



1 



.2a A \6= — 

2 r V roj 2m 2 



Q 2 



ezo* 



(21) 



(/ 2 -l) 2 + 



Q 2 



ezu* 



(22) 



(23) 
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We have defined the dimensionless parameters 



A77 2 



(24) 



where m\> is the mass of the gauge vector boson and Wd 
the mass of the dilaton in the Einstein frame for 



U(</>) = m^V/ 2 . 
The Ricci scalar is given by 



R 



-mi 4<r 



5a 2 



4(7- 



and w is the dimensionless angular metric factor 
vj = mhW. 



(25) 



(26) 



(27) 



The modulus of the Higgs field appears through the di- 
mensionless function f(p) 



V 



(28) 



Note that the metric (|15f) has no conical singularity in 
r = provided u> ~ r, i.e. tn ~ p in the hyperstring core. 
Differentiation of the action (J7J) with respect to the Higgs 
field yields the Klein-Gordon equation 



/ •-)./• & ~/(/ 2 !) 0. 

vj I vj 2 - 2 v ' 



(29) 



while differentiation with respect to the gauge field gives 
the Maxwell equation 



Q 



2& )Q-ef z Q 

w 



0. 



(30) 



Since the dilaton is only coupled to the metric, the 
Bianchi identities ensure that stress-energy conservation 
is still verified in the matter sector. Therefore, there is 
a redundant equation in the system of differential equa- 
tions (|20|l to H3()|l and we choose to solve numerically 
Eqs. (HU, and (gjty. 

As usual, the Higgs field vanishes in the defect core, 
i.e. $ = for p = 0, while it relaxes to its vacuum ex- 
pectation value (VEV), 77, in the bulk. These conditions 
translate into the following boundary conditions for the 
function /: 



/(0) - 0, 



lim / = 1. 



(31) 



The corresponding boundary conditions for the gauge 
field are given by 



0(0) = 1, 



lim Q = 0. 



(32) 



Moreover, we assume that no additional 5-like energy 
distribution lies in the hyperstring core. The geometry is 





FIG. 1: Transverse profiles of the dimensionless hyperstring 
forming Higgs / and gauge field Q, together with the metric 
factors e CT and wj p for the parameters a = 0.4, f3 = 0.1 and 
e = 0.1. The dilaton <j> condenses in the hyperstring core and 
the spacetime has an asymptotically flat (conical) geometry. 



therefore regular in the core and the metric coefficients 
verify 



a(0) 
vj(0) 



a(0) 
vj(0) 



0. 
1. 



(33) 



Asymptotically we are interested in the cosmic string 
branch solutions. Recall that the dilaton is massive and 
does not couple to the matter sector allowing the exis- 
tence of asymptotic flat conica l sp acetimes without long 
range dilatonic effects [M ES El El • This leads to the 
boundary conditions 



0(0) = 0, 



lim <f> = 0, 



(34) 



where the constraint at p = comes from the reg- 
ularity requirement of a non-winding field. The nu- 
merical method used to solve this system is described 
in appendix ^ and a typical solution is presented in 
Fig. n The dilaton condenses in the hyperstring core 
while the spacetime geometry is asymptotically conical, 
with a missing angle given by the asymptotic value of 
27r(l — tj/p) (see Fig. Note also the red/blue shift 
exp(cr) between the core and the outer regions. In the 
next section, the behaviour of the 4D tensor waves in this 
background is investigated, both analytically and numer- 
ically, and we will find that this system does not support 
quasi-localized gravity on the brane. 



B. Tensor Perturbations 

In this section, we investigate the propagation of ten- 
sor perturbations, with respect to the AD vortex world- 
volume, on the background field configuration studied 
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in the previous section. For this purpose we derive the 
equation of motion of tensor perturbations in the Jordan 
metric, assuming that the scalar and vector degrees of 
freedom are fixed. 

Since we are interested in pure AD tensor perturbations 
in the Jordan frame, the perturbed metric ljl5|) reduces 
to 

ds 2 = e° ( VliV + V) ^"dx" + aZ V Q " , (35) 
where /i M „ is a transverse and traceless tensor 

= °. t v K» = 0. ( 36 ) 
and z is a dimensionless conformal coordinate defined by 



-°^' 2 du. 



(37) 



With the hyperstring forming fields being of the scalar 
kind for the Higgs, and of the vector kind for the gauge 
boson, the tensor modes decouple from the background 
fields. The perturbed equations of motion are obtained 
by perturbing the Einstein-Jordan equations l|12l) at first 
order. From the perturbed metric tensor bg^y given in 
Eq. H35|) , one obtains in terms of the dimensionless quan- 
tities (see Appendix |B|) 

h "nu + [ l "' + — " + 0' ) C + + ^dghfiv = 0, 

\z vo J ^ 

(38) 

where the prime denotes the derivative with respect to 
z and □ stands for the dimensionless four-dimensional 
d'Alembertian 



- _ □ _ f^A 



(39) 



It is convenient to decompose these perturbations in 
terms of winding modes around the hyperstring and 
d'Alembertian eigenfunctions: 

M^.P I «)=E e *' / h<${M,p)A a (x*)dM, (40) 
p J 

where Ajg verifies 



(□ + M 2 ) A M = 0. 
Moreover, the rescaling 

AP) _ p 3<r/4+0/2 J— i (p) 



(41) 



(42) 



allows to recast Eq. (I38|l in a Schrodinger-like form for 

d 2 e 



the ^} modes 



dz 



a + V P (z)£ = M 2 C, 



(43) 



where the tensor and angular mode indexes have been 
removed. The potential V p is given by 



V p = W z 



w 



TP 



with 



4 I w I 



(44) 



(45) 



In Eq. Ij43(l . M 2 stands for the dimensionless 
d'Alembertian eigenvalue 



M 



M 2 



"T, 



(46) 



Note that for the zero angular momentum modes p = 0, 
Eq. 143 , |l is supersymmetric and W is the superpoten- 
tial |43,[44[- Indeed, defining the operator 



A= — + W(z), 



simplifies Eq. (|43|l into 



(47) 



(48) 



This ensures that all the eigenvalues M 2 are positive, and 
thus there are no tachyons in this model, provided super- 
symmetry is not broken, i.e. there exists a normalizable 
zero mode. 



C. Quasi-Localized States and 4D Gravity 

1. Pure Cylindrical Waves 

It is interesting to briefly recall the case of cylindrical 
waves propagating in flat spacetime. In that case, with- 
out dilaton, a = <f> = and w = z, the potential 144|) 
simplifies into 



p 



(49) 



which is always negative definite for p = and positive 
definite otherwise. Eq. 143|) is a Bessel equation whose 
regular solutions in the hyperstring core are 



C P >(M,z)oty/z~J p (Mz) 7 



(50) 



the Y p and K p modes being singular in z = whereas 
the I p modes are singular at infinity. Moreover, since |45j 



J p (z) 



o T(p 



— (-Y 



(51) 



only the zero angular momentum modes p — do not 
vanish in the string core. Note that in the special case 
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M = the regular solutions simplify to z p . Far from the 
hyperstring core 



2 / - 
— — cos Mz 
ttM \ 



(52) 



and using plane wave normalization at infinity, one gets 
from Eq. I|42(l the spectral density 



-M. 



(53) 



These tensor modes therefore contribute to a Newton po- 
tential on the brane [Til Il2| 



V s (x) cx 



h(M,0) 



-Mi 



dM oc 



(54) 



i.e. compatible with the standard Newtonian gravity in 
a six-dimensional flat spacetime. 

As discussed in the introduction, if there exist localized 
or quasi-localized waves in the hyperstring the spectral 
density will be modified, and hence so will the resulting 
potential in Eq. 



> -0.25 




FIG. 2: The potential z 2 V {z) governing the behaviour of 
the tensor perturbations around the hyperstring according to 
Eq. H4.'ifl ■ The background fields are those of FigQand the 
straight line represents the pure flat spacetime case. 



2. Gravitational Propagation on the Hyperstring 
Background 

According to the background solution computed in 
first section, a' = <f>' = and w' /m — 1/z both in the 
hyperstring core and at infinity (see Fig.^J. In these two 
regions the potential V p behaves like the flat spacetime 
one and only the zero angular momentum modes p = 
can contribute significantly to the spectral density. As 
can be seen in Fig. [21 significant deviations from the flat 
case appear in the intermediate region where the deriva- 
tives of the background fields are non- vanishing. 

However, the potential remains negative definite and 
does not allow for the existence of bound states or quasi- 
bound states in the hyperstring core. In that case, we 
can use the WKB method to approximate the solutions 
in the entire extra dimension since the turning points 
where M 2 = Vq do not exist (see Fig. |2J) 



exp 



£wKB ( M , Z ) « 



\J M 2 - V (u) du 



[M 2 -V (z)} 



1/4 



(55) 



After normalization at infinity one gets, from Eq. (|42|l . 
a spectral density in the hyperstring core similar to the 
one induced by the pure cylindrical waves 



Mo) 



(56) 



and thus a 6D Newtonian potential, albeit with a weaker 
gravitational coupling. 

We have plotted in Fig. [21 the values of h(M, 0) ob- 
tained from numerical integrations of Eq. (|43f) for a wide 



FIG. 3: Normalized values of h(M,0) as function of M (solid 
line). For high and small values of M, one recovers the \J~M 
behaviour (dashed line). As expected from the WKB approx- 
imation, only small deviations from the 6D gravity appear as 
long as the potential does not allow bound states. 



range of masses M. This confirms the WKB result, with 
only slight deviations in the intermediate range. 

According to the previous discussions, this is physically 
expected as long as the potential Vq remains negative def- 
inite: there are no metastable AD tensor waves inside the 
core that could induce peaks in the spectral density. In 
the next section, we search for the precise conditions un- 
der which the potential would in fact become positive for 
a general codimension N c = 2 defect and generic dilaton 
potential. 
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D. Confining Potentials and the Dominant Energy 
Condition 



According to Eqs. I(44|) and l(45|) . the potential Vo driv- 
ing the propagation of the non-winding (p = 0) AD tensor 
modes in the bulk reads 



1 /3 , 



zu 
zu 



2 ^'2 

w 



zu 



-41'. (57) 



From the topological defect stress tensor components 
(|19|l . one can define the dimensionless matter energy £ 
and pressure V by 



£ = 



1 



-4 



f 2 2 1 
3 zu 2 4 {J ' 



(58) 



The second order derivatives in Eq. I|57|l can be expressed 
in terms of first order derivatives by means of the dy- 
namical Einstein- Jordan equation (|20|l . In terms of the 
conformal coordinate z and U = exp (cr)[//m 2 , Eq. i|20|) 
yields 



t>" = -2ac a -*£ -lu- -a' 2 - a'— 



zu 



Making use of the constraint equation (|21|) 



2 2 



2ae" 



2a'- 



-P--U, 



2a' 



zu 
zu 



(59) 



(60) 



the potential Vq simplifies to 



2V Q = -2ae°-* {£-V)-\ (\a' + — 



-U, 



(61) 

which, according to Eq. 1)58(1 . is negative for positive dila- 
tonic potentials U . In fact it is clear in Eq. (|61[1 that, for 
a general defect in codimension two, either £ — V < 
or U < are required for Vq to become positive at some 
point. The former condition would violate the Domi- 
nant Energy Condition (DEC) for the matter stress ten- 
sor, and one would need "special" hyperstring forming 
fields that support tachyonic propagation 01- It is be- 
yond the scope of this work to discuss more fundamental 
motivations and consequences of negative dilatonic po- 
tentials however we note that the simplest exam- 
ple of this is achieved by choosing a negative cosmolog- 
ical constant [see Eq. J7J)]. In absence of a dilaton, the 
six-dimensional spacetime generated by a hyperstring in 
the presence of a negative cosmological constant has in 



fact already been studied in Refs. |29|, |48|, |49| and ex- 
plicitly realizes Randall- Sundrum gravity confinement. 
However, such defect models require strong fine-tunings 
on the model parameters, and stay in the realm of fi- 
nite volume extra-dimensions. It is worth pointing out 
that in this case, since the potential is asymptotically 
positive due to the negative cosmological constant, there 
is a normalizable zero mode (reminiscent of the trapped 
massless graviton of the Randall-Sundrum model) , which 
ensures that the potential Vo does not break supersym- 
metry. Interestingly, one may also expect to confine a 
discrete spectrum of massive gravitons. 

In the next section, we explore the conditions under 
which there exist quasi-localized 4D tensor waves on a 
classical defect in higher codimensions, concentrating on 
the case in which the extra dimensional space is isotropic. 



III. CLASSICAL DEFECTS IN ISOTROPIC 
EXTRA DIMENSIONS 

The action Q can be generalized in a straightforward 
way to a iV-dimensional spacetime and the Einstein- 
Jordan equations remain given by Eqs. (|T2|l and ltT5|l . 
The dynamics in the matter sector is driven by the La- 
grangian £ ma t which should involve the required fields 
and their interactions to allow the formation of a AD 
topological defect in N dimensions. In the simple case 
where the defect forming fields do not break the AD 
Poincare invariance. i.e. there is no current flowing along 
the defect |U HJ IH H3> thc AD components of the 
stress tensor verify 



(62) 



where £ is a function of the transverse coordinates only. 
Like the other components of the stress tensor, the pro- 
file of £ in the extra dimensions, and consequently the 
behaviour of the metric fluctuations 5g^ v , depends of the 
underlying model defined by £ m at- 

Here, we consider the generalization of E q. (1 151) to 
more than two isotropic extra dimensions [5J, l55| 

ds 2 = e^n^dx^dx" + e a(z) dz 2 + u 2 {z)dn 2 n , (63) 

where 



dn 2 n = J2 T ^<^ d02 > 



(64) 



is the interval of a n dimensional maximally symmetric 
space. Depending on the curvature k of this subspace, 
one may choose the coordinates such as Ti = 1 and 



) n sin2 

i=2 



J 3h 



(65) 



where T 2 is sin 2 (6*2), 2 or sinh 2 (6>2) for positive, null and 
vanishing curvature k, respectively. 
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In this case, the Einstein- Jordan equations l|12(l read 



-a +n — 

2 uj 



£ - -e a U 
2 

3 , 2 n{n — 1) uj 



e a n(n - 1) 



OJ 

' 2 3 



- 2 

12 



no | a 

UJ 



w' 



(66) 

for the (uv) components, while the constraint equation 
(zz) is 



1 



n(n — 1) uj 1 ' 



UJ 



,UJ 

2na' — 

UJ 



e a n(n - 1) 



2a' 



(67) 

with V the pressure along the z extra dimension. The 
tensor perturbations Sg^ — exp(cr)/i Alt/ around this met- 
ric end up being solutions of 



n— + (j)' ) ti + Uh^ y + ^hnh^ = 0, 



where L„ is a generalized "angular" differential operator, 
e.g. for k = 1 one recovers the hyperspherical Laplacian 



i— 1 



tan 6j 



(69) 



For the zero angular momentum modes, i.e. those which 
are non- vanishing on the brane for symmetry reasons, one 
recovers the Schrodinger equation (|43fl for the rescaled 
quantity 



(70) 



Using the above Einstein-Jordan equations, as described 
in Sect. Ill Dl the potential is found to be 



2V = -« 2 „e CT " 



2 I 2 



(8 - V) - e a U + -^n(n - l)k 

(71) 



UJ 



n- 



UJ 



In contrast to the hyperstring potential studied before 
[see Eq. (ffiTJl ]. it appears that for positive curvature k = 1 
the potential may now take positive values, whereas for 
k = 0,-1, one would still need to violate the DEC to 
allow for quasi-localized states .This result is not really 
surprising, as often in General Relativity, a standard way 
to mimic matter that violates positive energy conditions 
is to consider positive curvature space (e.g. see Refs. [2(], 
HE EE Ell)' From Eq. ll6^|) . note that k is the curvature 
of n — N c — 1 extra dimensional subspace and, according 
to the asymptotic behaviour of uj(z), the volume of the 
total N c extra dimensions may be infinite even for k — 1 . 



Although we have not presented the existence of soli- 
tonic solutions compatible with the metric H63JI , the con- 
dition k = 1 may be generically fulfilled. For example the 
case of N c — 3 extra dimensions with a positive curvature 
k may be realized by a hypermonopole [19L |5J, |59j . In 
such a framework one may indeed expect to find quasi- 
localized tensor modes in the core without violating the 

DEC m. 



IV. CONCLUSION 

We have discussed in this paper the conditions under 
which one can realize quasi-localized gravity using un- 
derlying topological defect models. For defects with codi- 
mension N c — 2 the existence of metastable tensor modes 
requires violation of the Dominant Energy Condition by 
the defect matter. For codimension iV c > 2 we found 
that this requirement may be relaxed for appropriate de- 
fect matter. In this case the background solution can, 
at least in principle, support metastable gravity. The ex- 
plicit construction of N c > 2 codimensional sigma models 
that quasi- localize gravity is subject to future work 60] . 
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APPENDIX A: RELAXATION METHOD 

To solve the set of differential equations ((201, 122(1 . 
(J2HJ), p9B and (|3()|l . we have used a relaxation 
method |61| . In terms of dimensionless fields and pa- 
rameters, the action (J7J) can be recast into 



2a 

7T?7 2 



(j) (47707 + 2wj 2 ) dp + 2a / 7 2 ?tj 



f 2 Q 2 



+ ,(/ 2 - 1 ) 2 + 



Q 2 



dp, 



(Al) 

with 7 = e CT , and where integrations by parts have been 
used to keep only first order derivatives in the metric 
factors. After discretization of the radial coordinates p, 
the discrete action reads S = S e + S m with 
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TTT] 2 ^r 1 



(ji+1/2 ~ Ti-l/a) , (li+l/2 - (™i+l/2 - ^i-l/2) , 2 (4+1/2 - ^i-l/2)' 
oTOj ; 47i h 7 4 



ft " ft ' ft 

(7i+l/2 _ 7i-l/2j (4+1/2 - 4-1/2) „ 2 ( ro «+l/2 _ ^i-l/z) (4+1/2 - 4-1/2J 



27f 



+ ft/?7i mifa 



in the gravity sector, and 



2 

7T7] 



E 



,2_ ■" 1 , hJj fjGi 1 u 7i ZU i tf2 i\2 , 7» 2 (Qj+1/2 ~ Qt-l/2) 

7i w i r 1" ^7* 1" n— r- (/* - 1) H r 

ft a7i 4 ezui ft 



(/i+1/2 ~ /i-1/2) 



(A2) 
(A3) 



for the matter sector. The grid resolution is given by ft — 
Pi+i ~ Pi where the index i indicates that the fields have 
to be evaluated at the discrete points pi. The discrete 
derivatives have been expressed in their centered form, 
e.g. 



fi 



+1/2 



fi- 



1/2 



0(ft 3 ), 



(A4) 



where /j+x/2 is evaluated on a half step shifted mesh. By 
differentiating the discrete action S with respect to the 
discrete fields 7$, Wi, (pi, fi and Qi, one gets the finite dif- 
ference equations corresponding to Eqs. (E2J), 
(|29|l and 130() . respectively. From an initial guess of all 
the discrete fields on the grid pi, the solutions of the 
finite difference equations are obtained by a successive 
over-relaxation method. Here Ti designs a five dimen- 
sional vector whose components are the discrete fields 
7i, zui, (pi, fi and Qi. At step p + 1, the fields are up- 
dated by a Newton's method to reduce the error with 
respect to the true solution: 

f! +1 =f? + s5ff, (A5) 

where SJ-f is computed to solve 

£ (Tf +1 ) ~ ) + W£(5fif) = 0, (A6) 

£ being the finite difference equations. In Eq. 1A5|I . s is 
the over-relaxation factor. The boundary conditions are 
part of the finite difference equations £ since they appear 
as the constraints which have to be satisfied by T\ and 
T V N , where N is the total number of points of the pi mesh. 
The iterative process is stopped when the discretized ac- 
tion S given by Eqs. I|A2(1 and (|A3|) remains stationary 
at the machine precision. Moreover, we have also verified 
that the relaxed solutions satisfy the constraint equation 
(I21|l . The relaxed fields have been plotted in Fig. \T\ior 
an assumed generic set of parameters a, /3 and e. 



APPENDIX B: PERTURBED QUANTITIES 

For the metric (|63|) , the only non- vanishing tensor per- 
turbations are 

Sg M u = e CT V, V = -c~ VVX/s, (Bl) 
from which the perturbed Christoffel symbols read 



8Y a 



1 



(B2) 



2 (K» + u> v) . 

1 e CT d ei h^ v 



ST"- =-- 

» u 2 to 2 T, 

From the perturbed Ricci tensor 

SR BD = -5T 



ST A +T H ST A -T A 6T H 

ui BD.A ~ A BD AH x DH U1 ~ BA 



1 BA DH • 1 Aff 01 BD1 



(B3) 

the non- vanishing components of the perturbed Einstein 
tensor read 



SG^ = -- Ih"^ + ( +n^- ) ft^ + Uh^ v 



2 L^/lyj,^ 



3cr" + 2n 



+n(n - 1)— j + 2na — 



oj 2 



n(n — l)fc 



(B4) 

From Eqs. Ijl2|l and (f B If) - the perturbed dilatonic source 
terms, i.e. other than the matter stress tensor, are 















('- 


fn- 











5^) 



(B5) 



10 



whereas 8T^ V is readily obtained from Eq. (|62|l . or and (|B5f) . one obtains Eas. I|38|) and l|68|l . up to the back- 
Eq. (|19|l for the Abelian Higgs vortex. From Eqs. (|U4[) . ground Einstein- Jordan equation 112|) . 
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